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The electronic structures of the actinide elements U, Np, Pu, Am, Cm and Bk are investigated 
within the self-interaction corrected local spin density approximation. This method allows to de- 
scribe a dual character of the 5/ electrons, some of which occupy localized and core-like states, 
while the remaining 5/ electrons hybridize and form bands. Based on energetics the calculations 
predict delocalization/paramagnetism in the early actinides, and localization/anti-ferromagnetism 
in the later actinides. The corresponding calculated equilibrium volumes are in agreement with the 
experimental values. For Pu and Am, the method wrongly predicts magnetic ordering, but we find 
that the paramagnetic state gives a better description of cohesive properties. Under compression, 
in the later actinides, a localization-delocalization transition happens gradually as more and more 
/ electrons become band-like with decreasing volume. Pu is already at this transition point at am- 
bient conditions. Delocalization sets in for Am and Bk at a compression of V ~ 0.75Vo, for Cm at 
V ~ O.6OV0, where Vb is the equilibrium volume, and the transition is complete for V ~ 0.4 — 0.5Vo 
in these three elements. 

PACS numbers: 71.15.Nc, 71.27.+a, 71.30.+h, 71. 23. An 



I. INTRODUCTION 

The electronic structure of actinides re- 
mains an active area of research both 
experimentallji^'^^'^'^'i'^'a^'ii^ and 
thcorctically .^'^-^^-^^-^^-^^-^^--^^-^"-^^'^^-^'^-^^-^^'^^-^''-^^-^^-^"-'^^ 
The intricate nature of the /-electrons, their sensitivity 
to external probes, like temperature and pressure, and 
their magnetic ordering are at the center of interest. 
Already in the seventies it emerged that with respect to 
volume behaviour one had to distinguish between the 
early actinide elements, with their bonding properties 
resembling those of the transition metals, and the late 
actinides, where the lattice parameters remain relatively 
unchanged, reminiscent of the localized /-states of the 
rare earth seriesi^ Pu is on the border line between these 
two opposing pictures with a particularly rich phase 
diagram.— These trends observed across the actinide 
series can be correlated with an increase of the on-site 
Coulomb correlation energy of the /-electrons which 
gradually conquers the band formation energy as one 
progresses through the series. With applied external 
pressure the late actinides exhibit volume collapses 
associated with structural phase transitions signaling 
the onset of /-bonding 

Most theoretical studies on the actinide 

j,^^,j^,,j^^ jl3.14.15.16.17.18.19.20.21.22.23.24.25.26.27.28.29.30.31.32 

are based on density functional theory (DFT), in either 
its local density or generalized gradient approximation, 
and their overall conclusion is that the early actinides, 
Th, Pa, U, and Np, are well described by the conven- 
tional band theory, as their /-electrons contribute to 
bonding, and the equilibrium structures and lattice 



parameters are in good agreement with experiments^^ 
However, for the later actinides the band picture fails, 
and additional assumptions or parameters derived 
from experiments need to be invoked, diminishing the 
predictive power of the approach. Thus, from Cm to 
Es, the lattice parameters are found to be in better 
agreement with experiment when the /-electrons are 
prevented from participating in bonding^ indicating 
the need of a different energy functional. For Pu 
and Am, situated at the localization-delocalization 
transition, neither of these two approaches is capable of 
describing their electronic and magnetic properties, in 
particular conventional band theory finds both Am and 
(5-Pu magnetic. The onset of spin-polarization in band 
theory may be interpreted as a signature of /-electron 
localizatioufi^ albeit only close to half-filling, and when 
combined with anti-ferromagnetic ordering and orbital 
polarizationiSiiiii^ a quite accurate account of bonding 
energetics is obtained. The predicted magnetic ordering 
is however in disagreement with the experimentally 
observed Pauli enhanced paramagnetic ground state in 
(5-Pu and Am. By allowing for a complete disorder of 
the local magnetic moments, Niklasson et alji^i obtained 
a reasonably accurate description of bond lengths and 
bulk moduli throughout the actinide series, however, the 
underlying assumption of fluctuating local moments has 
not been confirmed by experiments. 

The inadequacy of the local density functional for the 
later actinides can be traced to its inaccurate repre- 
sentation of strong electron-electron interaction, which 
tends to localize the /-electrons on their atomic sites. 
To overcome this shortcoming, all implementations of 
the LDA-I-J7 approach2£i^i2^iS^ introduce the effective 
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Coulomb U parameter that separates the /-manifold into 
lower and upper Hubbard bands, and removes /-degrees 
of freedom from the Fermi level. In the fully local- 
ized limit (FLL) LDA+U finds 6-Pu magnetic ^i^i2^ but 
a careful balancing of spin-orbit and exchange interac- 
tions in the implementation of Ref. [H leads to a non- 
magnetic ground state of 5-Pu. In the around mean-field 
(AMF) flavour of LDA-I-U non-magnetic ground states 
were found for both 5-Pu and Ami^i^ The recent de- 
velopment of dynamical mean-field theory (DMFT);^ 
combined with the LDA+U approach, has enriched 
the field in terms of phenomena accessible to calcula- 
tions, including both ground state cohesion^ phonons^^ 
and photoemissioui ^^i^^i^^ However, most applications to 
date invoke the Hubbard Hamiltonian and inherit the un- 
certainties associated with the LDA+U method, most no- 
tably the effective U and the double counting correction, 
shortcomings that might be overcome if DMFT could be 
successfully merged with the GW technique!^ The dis- 
ordered local moment approach of Ref. [l^ is related to 
DMFT by a static approximation. 

The self-interaction corrected local spin density (SIC- 
LSD) theory^! provides a dual picture of coexisting 
localized and band-like /-electrons, which was also 
conjectured^ to be the appropriate picture of the 6- 
phase of Pu. The method is fully ab-initio as both 
kinds of /-electrons are treated on an equal footing, with 
no adjustable parameters. Consequently, a unique to- 
tal energy functional can be applied to the entire range 
of actinide elements to determine the ground state lo- 
calized/delocalized /-electron manifold of each element. 
The method was previously applied to the actinide ele- 
ments Pu to Cf assuming a ferromagnetic ground state 
and, with the exception of Pu, the calculated volumes 
and bulk moduli were in reasonably good agreement with 
the experimental valuesi^ The motivation for the present 
paper is to extend the SIC-LSD study to other mag- 
netic orderings (para-, ferro- and antiferro-magnetism), 
and explore the consequences for the bonding proper- 
ties and valence state of the actinides. In addition, the 
delocalization of /-electrons under pressure in the later 
actinides is also investigated. Through the present more 
comprehensive study, we seek to establish and document 
the SIC-LSD applicability and status regarding a thor- 
ough physical understanding of such complex systems as 
actinides. 

The magnetic and paramagnetic states are realized 
by switching on or off the spin dependent part of the 
exchange-correlation potential. We find that, except for 
U and Np, the spin-polarized solution has always a lower 
energy than the non-magnetic solution, which in case of 
Pu and Am is in disagreement with experimental ob- 
servation. However, it emerges that assuming a param- 
agnetic representation of the localized states provides a 
much improved description of the equilibrium volumes of 
both 6-Pu and Am. Antifcrromagnetic order is found to 
be the lowest energy solution for both Cm and Bk, in 
agreement with experiment, providing slightly improved 



specific volumes compared to the values calculated for the 
ferromagnetic case. Since SIC-LSD describes well the de- 
localized nature of the /-electrons in the early actinides 
and the localized nature of the /-electrons in the late 
actinides, this demonstrates that the transition from /- 
bonding to /-localization in the 5/ series is reproduced, 
while an adequate description of the cross-over region 
depends strongly on the representation used for the lo- 
calized states. 

The paper is organized as follows. In Section II the 
important aspects of the SIC-LSD method are outlined. 
In Section III we present our results for the U, Np, Pu, 
Am, Cm and Bk metals. The observed trends of the 
present study are discussed in Section IV, and Section V 
gives our conclusions. 



II. THEORY 
A. The SIC-LSD total energy method 

The LSD approximation to exchange and correla- 
tion introduces an unphysical interaction of the electron 
with itself^ which, though insignificant for extended 
band states, may lead to uncontrollable errors in the 
description of atomic-like localized states, for example 
the /-electrons in the later actinides. The SIC-LSD 
method22i^i2S corrects for this self-interaction, by pos- 
tulating a manifold of coexisting localized and delocal- 
ized /-electrons, and by adding to the LSD total energy 
functional an explicit energy contribution for an electron 
to localize. The approach is fully ab-initio, as both lo- 
calized and delocalized states are expanded in the same 
set of basis functions, and are thus treated on an equal 
footing. 

One major advantage of the SIC-LSD energy func- 
tional is that it allows to realize and study different va- 
lence states of ions constituting a solid. By assuming 
atomic configurations with different total numbers of lo- 
calized states, self-consistent minimization of the total 
energy leads to different local minima of the same total 
energy functional, jtjSic-lsd ^ ^^^^ hence their total en- 
ergies may be compared. The configuration with the low- 
est energy defines the ground state configuration. If no 
localized states are assumed, E^^'^~^^^ coincides with 
the conventional LSD functional, i.e., the Kohn-Sham 
minimum of the E^^^ functional is also a local mini- 
mum of £:S/c-LSD Another advantage of the SIC-LSD 
scheme is the possibility to localize /-states of different 
symmetry. In the present work we exploit this possibility 
to investigate several magnetic structures including non- 
magnetic, ferromagnetic and antifcrromagnetic order. 

Given the total energy hmctional £'S'/c-Lsri^ -j-j^g 
computational procedure is as for the LSD case, i.e. 
minimization is accomplished by iteration until self- 
consistency. In the present work, the electron wavefunc- 
tions are expanded in linear- muffin-tin-orbital (LMTO) 
basis functions;^ and the energy minimization problem 
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becomes a non-linear optimization problem in the ex- 
pansion coefficients, which is only slightly more compli- 
cated for the SIC-LSD functional. Further details of the 
present implementation can be found in Ref. [stI . For 
a variety of a ppli cations to d and / electron solids, see 
Refs. l37ll4lH4al43ll44ll45l and references therein. 

The elementary excitations of solids, i.e. the one- 
particle energies, are often approximated by the band 
eigenenergies of the Kohn-Sham scheme, although a for- 
mal justification does not exist. Although for weakly 
correlated materials such an approximation may by rea- 
sonable, for localized electron states it is ill founded. It is 
possible to define band states constructed as Bloch com- 
binations of the localized states (sometimes referred to 
as 'canonical orbitals— ), however the associated band 
energies bear little resemblance to physical removal en- 
ergies, as e.g. observed in photoemission spectroscopy. 
This is in contrast to the LDA+U approach where bet- 
ter agreement with photoemission experiment can be ac- 
complished by an appropriate choice of the U parameter 
and double-counting correction. On the other hand, the 
SIC-LSD method, being fully ab-initio, does not require 
parameters, and can therefore be applied to study trends 
in energetics of the entire series of actinides in a com- 
prehensive manner, and with potentially predictive ca- 
pability. The SIC-LSD method is not a means to explain 
spectroscopies but should be treated as a useful total en- 
ergy approach, allowing to study a variety of phenomena 
and among them valences and valence transitions. 



B. Calculational Details 

Calculations have been performed for the face-centered 
cubic (fee), hexagonal close-packed (hep), and double 
hexagonal close-packed (dhcp) crystal structures. The 
Wannier states describing the localized / electrons are 
expanded on clusters of 32, 64 or 48 sites around the 
localization site^^ for these three cases. Similarly, the 
itinerant states have been sampled using 525, 216 or 75 
k-points in the irreducible wedge of the Brillouin zone, 
respectively. The 6s and 6p semicore states have been 
treated self-consistently in a separate energy panel. In 
the valence panel the LMTO basis set includes functions 
of 7s, 7p, 6d and 5/ character, where the 7p channel has 
been downfolded4S 

The atomic spheres approximation (ASA) is used, 
whereby the crystal volume is divided into slightly over- 
lapping atom-centered spheres of a total volume equal 
to the actual volume. The ASA gives rise to uncontrol- 
lable errors in the evaluation of the total energy, which 
inhibits the comparison of energetics of different crystal 
structures, on account of different overlaps of the atomic 
spheres. In the present study the focus is on the ener- 
getics of different localization scenarios within the same 
crystal structure, and the ASA error is of minor influence. 
The spin-orbit interaction couples the band Hamiltonian 
for the spin- up and spin-down channels, i.e. a doubled 



secular problem must be solved. Other relativistic ef- 
fects are automatically included by solving the scalar- 
relativistic radial equation inside spheres. The spin-orbit 
parameter, 

dr ' 

in atomic Rydberg units, is calculated from the self- 
consistent potential. 

Given the considerable disagreement between theory 
and experiment as to the magnetic ordering specifically 
in Pu and Am, we use the total energy considerations 
to investigate the SIC-LSD ground state with respect to 
both valency and spin configuration. The magnetic or- 
der is dictated by the order imposed on the localized 
/-states, and the paramagnetic and ferromagnetic/ anti- 
ferromagnetic scenarios are implemented by consider- 
ing the initial localized /-manifold to be either non- 
spin-polarized or spin-polarized in appropriate spatial 
arrangements.— During iteration to self-consistency the 
Wannier states may change, however roughly retaining 
their overall characteristics (exceptions occur close to de- 
localization) . 

Other approaches implement a non-magnetic ground 
state either by assumption (DMFT)^ through a care- 
ful calibration of the U and J parameters in LDA-I-U 
(FLL)^ or the double counting correction in LDA-I-U 
(AMF)f22 or by invoking disorder in the directions of lo- 
cal moments 



III. RESULTS 

As described in the previous section, we consider in 
this work only the high symmetry crystal structures rele- 
vant for the actinide elements in their low-pressure/large 
specific volume phases with localized /-electrons, specif- 
ically fee and hep crystal structures. A few tests for the 
dhcp structure, which is the actual ground state struc- 
ture at ambient conditions for Am, Cm and Bk, revealed 
virtually no difference with respect to the simpler hep 
structure in terms of cohesive energy per atom, specific 
volume and magnetic moments. As mentioned earlier, 
the present ASA implementation of the SIC-LSD energy 
functional is not accurate enough to resolve the energy 
differences between these structures. 



A. Total energy considerations 

In Table I, the results of our calculations are summa- 
rized. In the last column, the calculated total energy dif- 
ferences between the paramagnetic and magnetic config- 
urations are shown. We find that the non-magnetic con- 
figuration is energetically most favourable for the light 
actinides, U and Np, while all the remaining actinides 
prefer the spin-polarized ground state. Thus, from SIC- 
LSD total energy considerations we find the calculated 



4 



TABLE I: Equilibrium volumes, /-occupation and magnetic moments of U, Np, 5-Pu, Am, Cm, and Bk in the SIC-LSD 
ground state. Results for both the paramagnetic and the magnetic ground state are quoted. The numbers of localized and 
delocalized /-electrons are given, and for the moments both the spin and orbital moments (in Bohr magnetons, /is) are listed. 
All paramagnetic (PM) calculations were done for the fee structure, while ferromagnetic (FM) calculations were performed for 
the fee structure in the case of U and Np, and antiferromagnetic (AFM) calculations for the hep structure for Pu, Am, Cm, 
and Bk. Last column gives the energy difference (in eV) between the paramagnetic and magnetic ground states. 





Vb (a.u.) 


Ka^p(a.U.) 


Nf (loc) 


Af/(band) 


Ml 


2Ms 


AEpM-AFM (eV) 


u 


140.4 


138.9 





3.07 


0.0 


0.0 




Np 


122.4 


126.9 





4.30 


0.0 


0.0 




Pu (PM) 


163.8 


168 


4 


1.34 


0.0 


0.0 


1.14 


Pu (AFM) 


207.1 


5 


0.35 


-4.11 


5.25 


Am (PM) 


200.6 


198 


6 


0.40 


0.0 


0.0 


1.56 


Am (AFM) 


210.5 


6 


0.54 


-1.59 


6.45 


Cm (PM) 
Cm (AFM) 


171.0 
203.8 


202 


6 
7 


1.51 
0.26 




-0.15 



6.79 


2.29 


Bk (PM) 
Bk (AFM) 


205.1 
197.4 


189 


8 
8 


0.67 
0.43 



3.21 



5.36 


1.86 



ground states in agreement with experiment for U, Np, 
Cm, and Bk, but fail to reproduce the experimentally ob- 
served non-magnetic ground states of both Pu and Am. 

The fact that total energy considerations in LSD and 
also SIC-LSD favour the magnetic ground state in the 
actinide elements Pu and Am is a manifestation of the 
Hund's first rule rather than a consequence of interatomic 
exchange being significant. Neither LSD nor SIC-LSD 
can describe the appropriate multi-determinant eigen- 
states of L and S. This leads to an error in the en- 
ergy associated with local paramagnetic moments, since 
the exchange potential generated by the localized states 
and exerted on the non-/ electrons is determined by Ms 
(~ S). This is most clearly demonstrated for /®, where 
the correct ground state according to the Hund's rules 
in Russell-Saunders coupling is "^Fq, i.e. S ~ L = 3, 
J ~ 0. Hence, also Ms = 0, and accordingly no spin- 
dependent exchange-field is exerted on the conduction 
electrons by an actinide ion in the f^C^o) configuration. 
But LSD sets up a large spin-dependent potential, since it 
can only represent the ''^o ground state by Ms = 3. This 
description can be valid if the inter-atomic exchange in- 
teraction is strong enough for an ordered magnetic phase 
to be formed. However, it is an incorrect description of 
the cohesion of the conduction states, if the exchange 
interaction is not strong enough compared to the spin- 
orbit interaction which selects the J = ground state 
multiplet. The latter scenario occurs for Am. This also 
explains why the disordered local moment approach de- 
scribes the actinides so wellfi^ since the disorder (in the 
direction of spin quantization axis) leads to averaging of 
the spin-dependent part of the potential to zero. 

B. Uranium and Neptunium 

Studies of the electronic structure of U and Np, have 
demonstrated that the observed volumes and crystal 
structure can be explained by the conventional band 
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FIG. 1: (Color online) Total energy of fee Np as function 
of volume for several scenarios distinguished by the differ- 
ent number of localized /-electrons on Np ions. Both non- 
magnetic (spheres) and ferromagnetic (squares) configura- 
tions are considered. 



structure calculations (see for example Ref. [H and ref- 
erences therein). Nevertheless, despite this overall rather 
good agreement with experiment there are still discus- 
sions as to what band structure approach is best suited to 
describe volume trends f^^ii^i^ Furthermore, a compari- 
son between XPS data and a fully relativistic calculation 
has shown that correlations already have a measurable 
infiuence on the electronic structure of Npi^ 

In our calculations, we compare the band-/ scenario, 
where all / electrons are treated as itinerant, to a number 
of localized/delocalized /-electron scenarios, and find it 
to be energetically most favourable for both U and Np. 
In figure [1] the total energy versus volume is plotted for 



a number of Np /-electron configurations, both for the 
non-magnetic and the magnetic states. The global to- 
tal energy minimum is obtained for the f'^ configuration. 
Spin polarization is not energetically favourable, result- 
ing in a non-magnetic ground state in agreement with 
experimental data for both U and Np. This indicates 
that correlations are relatively weak in these elements, 
and confirms the adequacy of the LSD picture for treat- 
ing the electronic structure of the early actinides. Com- 
paring the calculated equilibrium volumes to the exper- 
imental values, we find very good agreement for both 
U and Np, but with some slight overbinding in Np, as 
is also observed in other band structure calculations. A 
comparison of the difference in total energies between the 
LSD ground state configuration and the nearest localized 
configuration gives respectively 0.38 cV for U and 0.20 
eV for Np, implying that the localized configuration is 
less unfavourable for Np, which indicates the increasing 
influence of correlations. 



C. Plutonium 

Prom the theory point of view, Pu is without doubt 
the most interesting element of the actinide series, as 
witnessed by the considerable number of papers already 
published on the subject. Its complex phase diagram, 
behaviour under pressure, and the absence of magnetism 
have still not been fully understood to date. Within DPT 
a number of theoretical studies have concentrated on the 
magnetic ground state, where spin-polarization mimics 
the effect of /-electron localization, which leads to cal- 
culated lattice parameters that are in rather good agree- 
ment with the experimental values, for both the a- and 
i5-phaseJi2iiS Other calculations have tried to take the 
experimentally observed non-magnetic ground state into 
account, either by assuming it as a starting point of their 
calculationiSi^ or by implementing it into their model by 
means of an adequate choice of parameters<^2i^i2^ Using 
the latter approach, LDA+C/ calculation o^^i^^ find a non- 
magnetic ground state for metallic Pu in both phases, 
driven by a strong spin-orbit coupling in the 5/-shell. 

In our calculations for Pu we merge elements of both 
of these approaches, as in what follows we assume a non- 
magnetic ground state, and compare to our previously 
published ferromagnetic calculations In figure [H the 
total energy of Pu in the fee phase is displayed for sev- 
eral localized paramagnetic /" configurations. Remark- 
ably, the energy minima for the scenarios with 0, 1, 2, 3 
or 4 localized states are nearly degenerate (within 0.03 
eV/atom), while localizing 5 /-electrons is marginally 
less favorable (0.11 eV/atom higher energy than the 
scenario) and localizing 6 /-electrons is 0.90 eV/atom 
higher in energy. The implication is that fee Pu is situ- 
ated virtually at the localization-delocalization transition 
at ambient pressure. The Pu /-manifold is dynamical, as 
the /-electrons may freely fluctuate between a localized 
/* shell and a fully itinerant state. It signals an elec- 
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FIG. 2: (Color online) Total energy of paramagnetic fee Pu 
as function of volume for several scenarios distinguished by 
the different number of localized /-electrons on Pu ions. 



tronic wavefunction of much more complex nature than 
can be realized with the SIC-LSD approach, which in the 
usual spirit of DPT seeks to describe the many-electron 
system with a single Slater determinant. The equilib- 
rium volume of the /^ state is calculated as 163.8 a.u., 
which is close to the experimental value for the S-Pu vol- 
ume of 168 a.u., while the calculated equilibrium volume 
of the /° state, 123.0 a.u., is close to, but smaller than, 
the experimental equilibrium volume of a-Pu, 135 a.u., 
quite a typical overbinding for the LSD approximation. 
Hence, it is suggestive that the /"* phase be stabilized by 
high temperature or alloying^ to form the (5-phase of Pu, 
while the low symmetry crystal structure of the a-phase 
most probably stabilizes the f'~' scenario. 

In figure [3] the density of states (DOS) of a number of 
these localization scenarios is displayed. In the LSD sce- 
nario of figure [3^, the /-manifold is split into a j = 5/2 
and a J = 7/2 band by the spin-orbit interaction. The 
exchange interaction is absent here, and the additional 
splitting within each of the two groups of bands is due 
to crystal field effects. Here the Fermi level is situated 
at the top of the filled / = 5/2 band, with a resulting 
large DOS at the Fermi level. In Fig. [S]], which repre- 
sent the DOS of the /^ paramagnetic configuration, four 
of the /-states treated as band states in Fig. have 
become localized, leaving the remaining /-states avail- 
able for band-formation. The width of the valence bands 
is 5 eV, however with most of the weight in the region 
from the Fermi level to 2.5 eV below the Fermi energy. 
A considerable hybridization of band /-states is seen in 
the range from -1 eV below to the Fermi level, including 
a high density of states at the Fermi level, in agreement 
with the experimental photoemission spectrum^ and the 
high specific heat^'* of 6-Pu. We note however that a 
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FIG. 3: Density of states of paramagnetic fee plutonium 
in a) the LSD scenario b) the loeahzed seenario and e) 
the loeahzed /® seenario. Full line is the total density of 
states, while dashed line gives the / partial density of states. 
Units are states per atom and per eV, and the energy is given 
relative to the Fermi level. 



photocmission spectrum is not directly comparable to a 
theoretical density of states, which strictly speaking only 
provides information on how the ground state is con- 
structed from one-particle states, while the final excited 
states play a crucial role in the photocmission spectrum. 



calculations, Refs. [20II21II23I also concluded on an lo- 
calized ground state. The recent LDA+U calculations 
of Ref. based on a reduced exchange interaction de- 
termined a non-magnetic ground state characterixed by 

ions, with a vanishing DOS at the Fermi level. A 
similar picture emerges from our SIC-LSD DOS for the 

configuration, shown in Fig. [SJ^. Energetically we 
find this configuration unfavourable by almost 1 eV, and 
the absence of any noticeable DOS at the Fermi level is 
in disagreement with experiment. The AMF flavour of 
LDA-I-U finds S-Pn in a mixed valent state^^ with a total 
/-occupancy of 5.44. In a recent DMFT calculation;^ 
Pu is similarly described as being in a mixed valence 
state (strong peak around the configuration), with an 
average /-occupation of 5.2, while another DMFT cal- 
culation using a different impurity solver suggest a value 
around 5.8.— The integrated number of band /-states in 
our /^ configuration is 1.34, therefore the total number 
of /-electrons in S-Pu comes out to be 5.34 (see Table 
I). This result is markedly different from the situation 
in the rare earths where the number of delocalized /- 



electrons never exceeds about 0.7 
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e. significantly 



more /-electrons are involved in the cohesion in S-Pn at 
ambient conditions. The total number of /-electrons is 
nearly a constant in the various localization scenarios /° 
- /5 (Nf = 5.48,5.42,5.36,5.36,5.34,5.39, respectively), 
i.e. what distinguishes these cases is the ratio of local- 
ized to delocalized /-electrons. The /^ scenario becomes 
energetically unfavorable (figure ^ precisely because it 
forces a high /-occupancy on Pu {Nf = 5.93, slightly 
below 6 due to the tails of localized states reaching into 
the neighbouring sites), which involves a relatively high 
intraatomic Coulomb interaction energy. 

The localized /-states are not shown in figure [3l since 
their SIC band energies are poor estimates of physical 
removal energies. A rough estimate of removal energies 
of the localized states may be obtained by the transition 
state argument,— which when applied to the /"* state of 
Pu places the localized states at -3.1 eV. In experiments^ 
there is no evidence of particularly strong emission at this 
energy, so this aspect of the electronic structure of S-Pu 
remains a puzzle. 

In conclusion, it emerges that the paramagnetic SIC- 
LSD implementation describes the fee phase of Pu sat- 
isfactorily, while the combination of the self-interaction 
correction and the exchange energy in the magnetic limit 
overestimates the tendency to localization. Experimen- 
tally, the vanishing of magnetic moments in Pu was re- 
cently re-confirmed Furthermore, electron energy loss 
spectroscopy branching ratio a^^i^^ clearly favor the jj- 
coupling picture over the LS-coupling (for a-Pu). 



In the earlier SIC-LSD calculation^, assuming a ferro- 
magnetic ordering, we found a Pu ground state comprised 
of localized /^ ions, with much too large a volume (30 % 
overestimate) , quite like the paramagnetic /^ scenario in 
figure [21 which has a minimum volume around 205 a.u.. 
By invoking ferromagnetic order in the LDA-I-U (FLL) 



D. Americium 

The SIC-LSD total energy curves for americium are 
displayed in figure 3) As for Pu we have investigated sev- 
eral localization scenarios from /*' to /*", assuming a para- 
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FIG. 4: (Color online) Total energy of paramagnetic fee Am 
as function of volume. Several paramagnetic configurations of 
the localized 5/ subshell are considered, as discussed in text. 
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FIG. 5: Density of states of hep americium in the param- 
agnetic state with 6 localized /-electrons per Am atom at 
V = Vo = 200.6 a.u.. Full line is the total density of states, 
while dashed line gives the / partial density of states. Units 
are as in figure 3. Only the itinerant electron states are shown. 
The localized states by the transition state argument may be 
estimated to be situated around -4.0 eV. 



magnetic phase. The lowest energy is found for the 
configuration witli an equilibrium volume of Vq ~ 200.6 
a.u. per atom, which matches perfectly with the experi- 
mental volume of Vexp = 198 a.u.. At this volume the 
configuration is 0.27 eV higher, and the Z*" configuration 
is 1.4 eV higher, i.e. the ground state appears quite 
robust at ambient conditions, contrary to the situation 
in Pu. This result is in agreement with recent photoe- 
mission studies^ that measure a solid ground state 



configuration. A non-magnetic ground state configu- 
ration has also been suggested on the basis of LDA+[/ 
(AMF) calculations ,H 

The density of band states of Am in the paramag- 
neitc ground state is displayed in Fig. O The cal- 
culation shown has been performed for the hep closed 
packed structure, but the density of states of the fee 
structure is very similar. The Am /-bands are almost 
exclusively of / = 7/2 character and situated around 1 
eV above the Fermi level, with a small tail extending be- 
low the Fermi level. The total number of /-electrons is 
Nf ~ 6.40, i.e. approximately 0.4 /-electrons reside in 
this tail. The localized states are not shown in the figure 
but the transition state estimate places them at an energy 
^ —4.0 eV. Compared to S-Pu, Am has approximately 
one less band / electron, slightly deeper localized states 
and more narrow unoccupied /-bands. Early photoemis- 
sion experiments^ on Am reveal the localized states as a 
plateau between 2 and 4 eV below the Fermi level, with 
traces of /^ characte r^^i^^ between -2 and -1 eV, how- 
ever with virtually no /-related emission at the Fermi 
level. Hence the electronic structure of Am is very differ- 
ent from that of (5-Pu. The photoemission spectra have 
been interpreted in terms of a mixed valent ground state 
of fluctuating /^ and /^ configurations dominated by 
the former. The present SIC-LSD ground state, com- 
prised of firm /^ ions coexisting with a tail of 0.4 band / 
electrons, can be considered a reasonable representation 
of the same mixed valent ground state. 

Upon compression the Am /-shell becomes more and 
more active in the bonding, as illustrated in figure 31 
where the energy curves of the configurations with fewer 
localized/more delocalized /-electrons come closer to the 
/^ curve. At a compression around V/Vq = 0.75 the 
/■^ curve becomes lower in energy, around V/Vq = 0.65 
the /^ curve dips below and becomes most favorable and 
so on, until around V/Vo = 0.52 eventually the /" curve 
becomes the lowest, meaning that the lowest energy is ob- 
tained with complete band-description of the /-electrons. 
Hence the /-delocalization happens gradually in the vol- 
ume range between V/Vq = 0.75 and V/Vq = 0.52. In 
terms of pressure, this corresponds to the range p = 16 
GPa to p = 35 GPa, as given by the common tangents of 
the respective energy curves. Of course, since the crystal 
structure in figure 2] is fixed to be fee this is not directly 
comparable to the experimental behavior of Am under 
pressure)^ where the /-electron delocalization is accom- 
panied by a sequence of structural changes. However, if 
we assume that the Ami (dhcp) and Amll (fee) phases 
are indicative of localized /-electrons, the experimental 
onset of delocalization is at the Amll— AmIII phase tran- 
sition, which occurs at V/Vq ~ 0.77 and p = 10.0 GPa. 
The /-electron delocalization is presumably nearly com- 
plete in the AmIV (primitive orthorhombic) , which is 
obtained at V/Vq = 0.63 and p = 17.5 GPa. This vol- 
ume is larger than the V/Vq = 0.52 deduced from the- 
ory, most probably due to the different crystal structure, 
as it is well known that /-electron delocalization favors 



8 



lower symmetry lattices^ and vice versa. The interpre- 
tation of the AmIII (face centered orthorhombic) as a 
structure of intermediate /-electron delocalization is also 
corroborated by total energy calculations r^, for which ad- 
ditional /-correlations (in the form of anti-ferromagnetic 
spin-polarization and orbital polarization) need to be in- 
corporated to achieve accurate description. Hence, there 
is qualitative agreement between the SIC-LSD theory of 
/ electron delocalization under pressure and experimen- 
tal observations for americium. 

In our earlier SIC-LSD study of amcriciumSi only the 
ferromagnetic spin arrangement was considered. The cal- 
culated ground state configuration was also /^, however 
with an 8 % too large equilibrium volume, and the f"^ 
configuration only 0.12 eV higher in energy. In addi- 
tion, the spin and orbital moments were far from can- 
celling {Ms = SAlfiB, Ml = — 1.58/is, wher is a 
Bohr magneton). The situation is slightly improved by 
allowing for an antifcrromagnctic spin arrangement. In 
the hep structure, the energy per atom is lowered by 
0.03 eV in the antifcrromagnctic structure compared to 
the ferromagnetic arrangement, and the equilibrium vol- 
ume is lowered by 5 %, but the moments remain almost 
unchanged {Ms = 3.23^b, Ml = — 1.59^b). Hence, 
antifcrromagnctic ordering constitutes the lowest energy 
solution of the SIC-LSD minimization. Overall, however, 
the paramagnetic scenario provides the best description 
of the elemental americium, while imposing an antifcrro- 
magnctic order can only be considered as a simple means 
of partially patching up the bonding properties. A SIC- 
LSD study on Am pnictides and chalcogenides predicted 
large DOS at the Fermi level, in agreement with suscep- 
tibility measurements^, but not observed in photoemis- 
sion studiesi^ The results were obtained based on the 
assumption of a magnetic ordering, but invoking instead 
a paramagnetic state might provide some new insight into 
their electronic structure. 



E. Curium 

The SIC-LSD total energy of curium is shown for both 
the magnetic and paramagnetic scenarios in figure [B] For 
simplicity, only the fee structure is considered, except 
for a few calculations performed with antiferromagnetic 
spin arrangement in the hep structure (asterisks in fig- 
ure [6]). The half- filled /^ configuration of localized states 
is seen to give the minimum, at a volume V — 203.8 
a.u., which is in good agreement with the experimen- 
tal equilibrium volume. From Table I, and figure [SI we 
find the paramagnetic scenario unfavourable in energy 
with respect to the anti-ferromagnetic scenario by 2.3 
eV. Here the magnetic state with localized /^ ions hav- 
ing S = 7 /2, L = 0, and J = 7/2, provides an adequate 
description of Cm. The paramagnetic state in this case 
gives a poorer description, with /^ localized / shells and 
a 17 % too small volume. The antifcrromagnctic state is 
lower than the ferromagnetic by 0.04 cV per atom in the 
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FIG. 6: (Color online) Total energy of fee Cm as funetion of 
volume for several seenarios differing with respeet to the num- 
ber of loealized /-eleetrons on Cm, notably 6, 7, or 8 localized 
/-electrons. The balls correspond to the ferromagnetic order- 
ing scenario, while triangles correspond to the non-magnetic 
scenario, as discussed in text. In addition, the asterisks cor- 
respond to hep antiferromagnetic arrangement of /^ curium 
atoms. 



hep structure, but we are not able to resolve the struc- 
tural energy difference between the very similar fee, hep 
and dhcp phases. The antiferromagnetic phase is also 
observed expcrimentallyj^ A recent DMFT— study has 
similarly established the /'' atomic configuration for Cm. 

The relative stability of the half-filled /^ shell is ev- 
idenced by the high compression {V/Vq ~ 0.60) needed 
for the next delocalization step, /^, to be of comparable 
energy. This is also in accord with the higher pressure 
needed experimentally to convert Cm into a low sym- 
metry structure indicative of /-contribution to bonding 
(Cmlll monoclinic structure), appearing at p = 37 GPa 
and V/Vo = 0.65^ 



F. Berkelium 

The SIC-LSD total energies as a function of volume, 
calculated for berkelium in the fee structure and the fer- 
romagnetic spin- alignment, are shown in figure [T) The 
localization /" scenarios, with n = 0,1, ..,8, have been 
considered, for both magnetic and paramagnetic states, 
but only the magnetic configurations are shown in the fig- 
ure. The ground state is found for an /* localized shell, 
at a volume V ~ 201 A a.u., which is 6 % larger than 
the experimental equilibrium volume. The energy dif- 
ference between the ferromagnetic and antiferromagnetic 
arrangements in the hep structure is quite small, ~ 0.02 
eV per atom, in favor of the latter. The equilibrium vol- 
ume at the same time decreases to 197.4 a.u., which is a 
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FIG. 7: (Color online) Total energy of ferromagnetic fee Bk as 
function of volume for several scenarios differing with respect 
to the number of localized /-electrons on Bk. 
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FIG. 8: (Color online) Trends in the Wigner-Seitz radius 
through the actinide series. The (red) circles and (black) 
squares are calculated (in fee structure) and experimental 
Wigner-Seitz radii, respectively. For Pu the f* non-magnetic 
ground state is taken and compared with the experimental 
Wigner-Seitz radius of (S-Pu. The calculated (/") and exper- 
imental a-Pu Wigner-Seitz radii are marked with an open 
circle and an open square, respectively. 



small improvement over the ferromagnetic volume. The 
number of band /-electrons is 0.43. The magnetic mo- 
ments of the Bk atom are large, Ml = 3.21^b, 2Ms = 
b.iQfiB, leading to a total moment of M = 8.57/is. The 
experimental moment of Bk is even higher, ^gg? = 9.8/iB 
(Ref. [H), leading to a projected value M ~ 9.05, i.e. 
very close to that of a pure p Hund's rule coupled ion 
(here we assumed a ^-factor of 1.5 as appropriate for p). 

The configuration with an /*" localized shell is 0.23 
eV higher in energy. This state has 1.45 band-like /- 
electrons. Upon compression this scenario becomes pro- 
gressively more important, and at V/Vq — 0.75 it be- 
comes lower in energy. The high stability of the /^ shell is 
manifested through this remaining the lowest energy sce- 
nario up to compressions of V/Vq ~ 0.40, after which the 
/ manifold quickly delocalizes. The transition pressure, 
corresponding to the common tangent between the /^ 
and p total energy curves, is p = 15 GPa. Experimen- 
tally, Bk is observed in the dhcp phase under pressures 
up to p = 7 — 8 GPa, after which the fee phase occurs and 
is stable up to p = 25 GPa and V/Vq = 0.70,^ At pres- 
sures above 25 GPa a low symmetry phase, possibly of 
a-U type, is observed. The conventional interpretation^ 
is that the first two highly symmetric crystal structures 
are associated with localized /-electrons, and hence that 
/-delocalization sets in at 25 GPa, which is somewhat 
higher than the theoretical value found here, but accept- 
able within our limits of accuracy. It would be interesting 
to investigate whether signatures of a state characterized 
by /-electrons in a mixed /*" localized and 1.5 / itiner- 
ant picture could be observed for berkelium. Such a state 
would have similarity to the (S-Pu phase, i.e. both distinct 
localized features, like deep-lying /-states in photoemis- 
sion, and itinerant characteristics, like a large density of 



states at the Fermi level. However, in contrast to (5-Pu 
the Bk atoms would possess large magnetic moments, 
arranged antiferromagnetically. The present study is re- 
stricted to the high-symmetry structures fee, hep and 
dhcp, and thus it cannot be excluded that around the 
/-delocalization transition the energy gain by full delo- 
calization and adoption of an appropriate low-symmetry 
phase outweighs that of only partial delocalization. 

The paramagnetic scenario was also investigated for 
Bk (not shown). Energetically, this state is unfavorable 
(see Table I) , but also other aspects of the calculation are 
in disagreement with experiment. The /^ localized shell 
is found to be the lowest energy state in this case, with 
an equilibrium volume oiV = 205.1 a.u., i.e. somewhat 
larger than the volume of the magnetic phase. However 
the /^ localized scenario is very close in energy (0.02 
eV per atom), implying a possible low-pressure (p ~ 1 
GPa) transition to this state, which furthermore would 
be accompanied by a ~ 10 % volume discontinuity. Ex- 
perimentally, there is clear evidence of large antiferro- 
magnetically ordered moments 

IV. DISCUSSION 

Table I and figures [U [U and [10] summarize the find- 
ings of the present paper. The rather good agreement 
between calculated and experimental Wigner-Seitz radii 
is demonstrated in Fig. [8] Figure [9] shows the calculated 
number of band-like /-electrons in the ground state. It is 
seen to rise steeply from U to a-Pu, reflecting the filling 
of the /-band, before dropping dramatically towards the 
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FIG. 9: Trends in the the number of band / electrons through 
the actinide series. For Pu, the a- and (5-phases are identified 
with the /" and paramagnetic states, respectively. 



heavier actinides, reaching its lowest value for Cm, after 
which a small increase is seen again for Bk. For (5-Pu, an 
intermediate /-band filling is found, which elucidates the 
extraordinary character of this element, also clear from 
the essentially vanishing localization/dclocalization en- 
ergy illustrated in figure [2l Correlating Fig. [8] to Fig. [9l 
it follows that the increasing number of band-like / elec- 
trons for U, Np and Pu(/°) contributes to the bonding 
and leads to the observed smaller equilibrium volumes, 
while larger volumes are obtained for Am, Cm, and Bk, 
having a few band-like / electrons. 

Several important energy scales may be identified for 
the actinide elements. The intraatomic exchange field 
is undoubtedly relevant but cannot be comprehensively 
addressed by the single-determinant description underly- 
ing the DFT, by which the atomic S and Ms quantum 
numbers get mixed up. In our calculations this shortcom- 
ing translates into wrong predictions of magnetic ground 
state configurations of Pu and Am. The energy scale for 
the interatomic exchange interaction is small, favoring 
the antiferromagnetic arrangement over the ferromag- 
netic arrangement with 0.02-0.04 eV per atom in Pu, 
Am, Cm and Bk. The localization/dclocalization energy 
scale is larger but depends strongly on volumci^ As a 
consequence, /-electrons are eventually transferred from 
being localized to becoming band-like when pressure is 
applied. In figure [TOl we show the trends in the range of 
relative volume over which the /-electron delocalization 
takes place for Pu to Bk. The volume range where the 
/-electron delocalization occurs is defined in theory by 
the lowest volume having the same number of localized 
electrons as in the ground state (onset), and the largest 
volume where a completely delocalized (/'^) /-manifold 
occurs. In experiment there is no clear-cut definition, but 
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FIG. 10: (Color online) Trends in the ranges of volumes of / 
electron delocalization in the actinide elements. Onset (com- 
pletion) of delocalization is marked with balls (squares), for 
each element with experimental data (blue and full line) to 
the left and theoretical data (red and dashed line) to the right. 
Volumes are given relative to the (localized) equilibrium vol- 
ume at ambient conditions. Experimental ranges are defined 
by the smallest volumes observed in the high-symmetry (fee) 
phase and the largest volumes observed for the low-symmetry 
(a-U type or similar) phases in high pressure experiments 
(Refs. Ilo, [Til . 113; for Pu the range is defined by the zero pres- 
sure volumes of the 5- and a-phases) . The theoretical ranges 
are defined within the fee structure only. 



presumably as long as the crystal phase remains fee under 
compression, the /-shell is completely localized, whereas 
a complete delocalization may be assumed to have oc- 
curred once the high-pressure low-symmetry phases of 
AmlV^ CmlV^ and Bklll^ are reached. The trends 
of experiments and theory are in agreement, with Pu at 
the delocalization threshold under ambient conditions, 
and Cm the most stable of the actinide elements studied. 
The theoretical volumes of / delocalization completion 
tend to be smaller than the experimental ones, which 
most probably is due to the restriction of the present cal- 
culations to the high-symmetry fee and hexagonal phases. 

A striking result of the calculations presented here is 
that an excellent account of the cohesive properties of 
the localized phases of the actinide elements can only be 
achieved, if S-Pu and Am are described by a paramag- 
netic phase, while a magnetic phase is appropriate for 
Cm and Bk. The SIC-LSD total energy functional incor- 
rectly favors the formation of a magnetic state for (5-Pu 
and Am. The reason for this shortcoming is probably 
different for these two elements. In Pu, the /-shell is in a 
vividly fluctuating state, and presumably any tendency 
to form interatomic magnetic correlations is washed away 
by quantum fluctuations. In Am, the /^ shell is more 
robust against fluctuations, but here the overestimated 
tendency to form a magnetic state can be traced back to 
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the inadequacy of the LSD to describe Russell-Saunders 
coupling and the formation of Hund's rules ground states, 
as discussed in section III. A. 

Due to its construction, the SIC-LSD can only describe 
an integer number of localized states on each atom, and 
the delocalization process upon compression is modelled 
as a series of steps, characterized by one more /-electron 
changing from atomic to band-like behavior. Such a de- 
scription can only be considered approximative, and in a 
more elaborate theory the transitions are likely to occur 
continuously. To describe these continuous transitions 
at finite temperature we can construct a mean-field co- 
herent potential approximation of the fluctuating local- 
ization configurations distinguished by different numbers 
of localized states. This has already been implemented 
for cerium^ and a similar study for the actinides is in 
progress. Bk may be an exception, due to the high sta- 
bility of the shell, and one may envisage a pressure 
range where the shell is still localized while the eighth 
/ electron, of spin minority character, has delocalized. 
This would be an interesting effect to look for. A more 
complete study of the delocalization process would also 
need to consider the effects of crystal structure, as it is 
well known that /-electrons once starting to form bands 
also prefer low-symmetry structures. This necessitates a 
full-potential description of the total energy, which will 
be subject of future development of the SIC-LSD method. 

Finally, the most important point of Fig. 10, namely 
Pu situated at the delocalisation threshold under ambi- 
ent conditions, and trivalent Cm to be the most stable 
actinide studied here, is confirmed by a recent DMFT 
calculation^ This SIC-LSD result shows that the SIC 
approach is capable of identifying the important valence 
states across which the DMFT has to fluctuate dynam- 
ically: more or less all possible valence fluctuations for 
Pu, and the pure trivalent state for Cm. 

V. SUMMARY 

The cohesive and electronic properties of the ac- 
tinide elements U, Np, Pu, Am, Cm and Bk have 
been investigated with the self-interaction corrected lo- 
cal spin density approximations. The predicted non- 



magnetic/magnetic ground states are in agreement with 
experiment for the light actinides U, and Np, and the 
later actinides Cm and Bk, but a magnetic ground state 
is wrongly predicted for Am and Pu. However, an ac- 
curate description has been obtained (see Table I) , when 
Pu and Am are forced into a paramagnetic state. Cm 
and Bk are correctly described in an antiferromagnetic 
arrangement. U and Np prefer a non-magnetic delocal- 
ized / ground state configuration. The ground states at 
zero temperature and pressure are the trivalent Am(/^), 
Cm(/'^), and Bk(/^), while Pu shows almost degener- 
acy between any of the localization configurations from 
/° to /^. The latter signals a true ground state of a 
much more complex nature than can be described with 
the present approach, characterized by quantum fluctu- 
ations between any of these states. If the stabilization 
of Pu(/'*) with high temperature or Ga alloying is en- 
visaged, this state may be identified with the (5-phase of 
Pu, with an intermediate number of localized and itin- 
erant /-electrons and an accompanying large density of 
states at the Fermi level. The delocalization of the /-shell 
under compression has been studied, and the delocaliza- 
tion has been found to occur over a range of volumes, 
with the established trends in good agreement with ex- 
periments. However, to describe the actual high pressure 
crystal sructures and the structural transitions, a more 
accurate implementation of SIC-LSD, most notably in- 
cluding full non- spherical potentials, is necessary. 
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